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1 Introduction 



In ref. |l|] we studied the problem of transmission and reflection of fermions through 
a wall and established the relationship which connects the complex plane poles of 
the wall profile function with the high energy reflection coefficient. The aim of the 
present paper is to study in more details this previous result and to go further into its 
consequences, mainly CP violation. The interest in the above-mentioned old problem 
has been renewed in the last years by the proposal that the baryon asymmetry of 
the universe might have been produced if a first order SU{2) x U{1) phase transition 
took place during the cosmological evolution. This electroweak phase transition can 
be described in terms of bubbles of "true" vacuum with an inner expectation value 
of the Higgs field v ^ 0, i.e. a spontaneously non-symmetric phase, appearing and 
expanding in the preexisting "false" vacuum with v = 0. If the weak phase transition 
is first order the standard model of weak interactions has qualitatively all the ingredi- 
ents to produce the cosmological baryon asymmetry 0, i.e. baryon number violation, 
C and CP- violation and out-of-equilibrium processes 0]. The question is whether it 
can produce a large enough asymmetry to account for the observed baryon num- 
ber of the universe: ub/s ~ (4 — 6)10~^^. CP-violation in the Kobayashi-Maskawa 
scheme 0] has been claimed to be large enough to generate alone the observed baryon 
asymmetry . However a zero temperature estimate [^] of the GIM suppression of 
electroweak C and CP effects, as well as a finite temperature study stressing the effect 
of gluonic thermal fluctuation [0 , allows to argue that the SM scenario |^ produces 
a CP asymmetry more than ten order of magnitudes too small, ?t-b/s<10^^'^. 

Simple extensions of the SM may produce much larger CP asymmetries, such 
as the two-Higgs-doublet [§|, 0, |10[ and |Tl|, left-right symmetric models [jl2|, SUSY 



T3| or models with heavy leptons fl^, [jT5|. These models contain several Higgs flelds 



and a relative phase between the Higgs vacuum expectation values may generate a 
new source of CP violation. The latter CP violation may be much larger than in 
the standard model for two reasons: i) CP violation often appears at the tree level 
while in SM it is a loop effect; ii) there is more freedom in the parameters than in the 
SM. Even so, the parameters in the model of ref. ||15| are too constrained to allow 
for enough CP violation, but the discrepancy is only two orders of magnitudes as 
compared to the ten orders of magnitude of the SM. 

We will therefore concentrate on tree level CP violation due to vacuum ex- 
pectation values of Higgs flelds beyond the standard model. These CP-violating 
processes do not directly violate baryon number nor fermion number. But they can 
create local inhomogeneities in the baryon (lepton) number, for example when the 
fermions in the plasma hit a bubble wall: an excess of antibaryons (baryons) may be 
created outside (inside) the bubble, which may subsequently be converted [0, via 
anomalous weak interaction fl^, into a baryon asymmetry in the universe. 

We thus need to study the problem of particles, propagating in the plasma. 



say from outside the bubble, penetrating the wall and being eventually reflected out- 
side or transmitted inside the bubble. This is a difficult problem since the fermions 
interact not only with the wall (i.e. the Higgs vacuum expectation value) but also 
with the particles in the surrounding plasma. It involves solving a quantum Fokker- 
Planck equation taking into account CP-violation and baryon anomaly. A useful 
simplifying assumption has been proposed, the so-called "charge transport" mecha- 
nism PI, in which the anomalous weak interaction happens at a distance from the 
bubble wall. This assumption consists in decomposing the process into two steps, 
one describing the production of the CP asymmetry when the quarks/antiquarks are 
reflected/transmitted on the wall, the second describing the transport and eventual 
transformation of this CP asymmetry via the baryon number anomaly. We will con- 
centrate on the first step, the interaction of the incoming fermion with the wall. It 
may be assumed that during this first step the diffusion corrections are relatively 
minor ones to the scattering from the bubble wall, which is driven by the Higgs ex- 
pectation value (during the second step they are, of course, crucial). The thermal 
effects are only considered via the introduction of the Higgs field effective potential, 
which takes into account the thermal bath temperature. For simplicity we will thus 
describe the reflection/transmission mechanism of fermions at zero temperature on a 
wall generated by a Higgs expectation value which depends on x. The thermal effects 
could also be present via the thermal modification of the fermion spectrum (thermal 
masses). We will not consider them, although we believe that our conclusions in the 
present paper can easily be extended to incorporate them. Furthermore, the infor- 
mation given for the solution of the academic T = case can also be of use in other 
problems of physics, in solid state physics, in astrophysics, etc. 

The study of the reflection and transmission of fermions on a bubble surface 
is thus crucial to elucidate whether it may produce a large enough CP-asymmetry 
to generate the correct baryon asymmetry of the universe. The size of the bubble is 
in general very large compared to the typical size of the CP violating process under 
study. Therefore we will approximate the bubble surface by a fiat wall perpendicular 
to the z direction. The Higgs expectation values depend on z, and the resulting 
effective mass of the fermion is a function of z. This function depends on the above 
mentioned Higgs field effective potential. The profile obtained by solving the equation 
of motion is therefore rather complex and depends on many coupling constants [Q . 

A frequently used simplifying assumption is to represent the wall profile by a 
step function considering an extremely sharp phase transition which allows to com- 
pute, in an exact way, the two-point Green function for a free fermion in the presence 
of this thin wall ||^. However, it is known and will clearly appear later in this paper 
that in such a thin wall approximation all CP asymmetries due to Higgs expectation 
value at the tree level do vanish. The CP violating phase can just be rotated away 
by a chiral rotation. Therefore we will now concentrate on thicker walls, although 
the wall thickness must be smaller enough than the mean free path of the fermions 
in the plasma to neglect the diffusion corrections. 



In general, for more complex wall profiles the task of obtaining the two-point 
Green function is very difficult and analytic results for this Green function may not 
exist in most cases or be in practice unreachable. One may try numerical calculations, 
but, since Green functions are not really necessary at tree level we will not insist 
in this direction. Still, the first step is to find an orthonormal and complete set 
of eigenf unctions. The orthogonality of these eigenf unctions is a crucial question 
which can be considered in a general way. In appendix B, we show the orthogonality 
of the two different types of eigenstates built by requiring appropriate asymptotic 
behaviours far from the wall as in refs. |jl8|, ||l9l. One describes the wall by an ansatz 



that simulates the dynamics of the phase transition |Q [|T8[ [|T9[. There exists one 
such ansatz for which, in refs. [|18| and |19[, an analytical solution is given and a 



numerical calculation has been performed in |Q]. In both cases a new net effect of 
CP-violation is obtained at tree level, that cannot be obtained for the thin wall. 

It is known that observable CP asymmetries can only be obtained if, besides 
a CP violating (CP-odd) phase, a CP-even phase (a phase which is equal for parti- 
cles and anti-particles) exists and interferes with the CP-odd one. In the thin wall 
approximation a CP-even phase only appears in the total refiection domain. This is 
no longer true in the thick wall case, and CP-violating effects may exist outside the 
total reflection region, up to inflnite energy of the incoming fermion. 

A completely general solution of the thick wall problem is difficult. However, 
there is a domain in which a few exact results can be derived: when the energy of 
the incoming fermion is large. In the present paper we develop a general method 
to calculate in the high energy limit the refiection coefficient of fermions hitting a 
wall, establishing a simple relationship between the analytic properties of the profile 
function and the high energy behaviour. The quantum correction to the expected 
classical behaviour is obtained and its importance depending on the above mentioned 
analytic properties is showed. As we will point out, the range of energy for the quarks 
in the electroweak phase transition is appropriate to consider this limit, except for the 
top. The CP-violating effects are incorporated by means of an appropriate imaginary 
mass term and the basic quantity \R{E)\'^ — \R{E)\'^ is obtained. Interesting general 
theoretical consequences for the fermions quantum scattering and for the CP-violating 
effects will be derived from this connection of the analytic properties of the profile 
function and the high energy behaviour. Apart of the clear interest of this problem 
in cosmology, it is possible to apply the formalism to systems of relativistic fermions 
that can suffer a phase transition. Two examples can be condensed matter under 
extreme external conditions or certain stages of the quark-gluon plasma formation 
process. 

In section 2 we formulate the problem and obtain an integral expansion in the 
parameter rrio/E, where mo is the height of the wall. We consider separately analytic 
and non-analytic functions as wall profiles. Section 3 is devoted to the study of 
functions defined as zero outside a certain region, in which the wall is located, where 
the potential varies. These functions must be non-analytic on the real axis as will be 



seen. In section 4 general analytic functions are treated. An imaginary, chirally odd 
and CP violating mass term is incorporated into the formalism in section 5. Finally, 
in section 6 we summarize and conclude. 



2 The integral expansion for the Dirac equation 
solutions 



As usual in this kind of problems, we work in the rest frame of the wall, which is 
taken to be planar and parallel to the x-y plane. Thus, the mass will be character- 
ized by a z- dependent function, the wall profile function. The approximation for a 
planar interface for the bubble wall should be valid for large bubbles compared to 
the microscopic size scale. In the electroweak phase transition this is valid for most 
of the evolution of bubbles pOj. For simplicity we study the Dirac equation for one 
flavour. In order to calculate the reflection coefficient we need only the plane wave 
solution for particles moving along the z-axis. In any case, general solutions for other 
incoming directions can be obtained by performing the appropriate Lorentz boost in 
the x-y plane. 

Following Nelson et al. 



141, we work in the chiral basis and factor the Dirac 



operator into 2x2 blocks. Thus the Dirac equation can be expressed as 



(id, + idt -m*{z) \ 

m{z) idz — idt 

idz + idt —m{z] 

V m*{z) id, -idt J 
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(1) 



and using the following ansatz for solutions with positive energy E 
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we obtain 
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Taking into account that the eigenvectors of the chiral basis have been reordered to 

.5=f^3 

-0-3. 



obtain 7^ 







the eigenvalues of this matrix, the chirality, will be + 1 for 



ipi and '?/'4, and - 1 for ip2 and ip^. 



The solutions of the equation (^ can be written as follows 



(5) 



Where V indicates a path ordered product and r is the position variable along the 
z-axis. The Dirac equation solutions are determined by 



n{z,zo)=Ve'^^o''^'^^^^^ and 



which, in a matrix way, can be written as 



n{z, zo) 



UJ21 UJ22 



Q{z,zo] 



UJ21 UJ22 



(6) 



(7) 



It is very easy to show that 0721 



LU 



12 



and UJ22 



LU 



11 



and analogously for Q{z,zq). 



Thus, it follows from (H) that, knowing Un, UJ12, cJn and cJi2, we have the solutions 
of the time-independent Dirac equation . 

If we consider 7n{T) = for t < Zq, then ipi and 1^2 correspond to right-moving 
right-handed particles and left-moving left-handed particles respectively, as well as 
ips and ip4 to right-moving left-handed and left-moving right-handed particles, in 
this region. If m(r) is considered to be a certain constant, m, for t > z, r2(r, z) and 
Jl(r, z) can be immediately diagonalized in order to identify the right-moving and the 
left-moving flux of particles. Thus, by requiring that the left-moving flux is zero for 
T > z, i.e. allowing only a transmitted flux beyond the wall, a left-handed reflected 
flux is obtained from the right-handed incident one and vice versa. Therefore 



n 



R^L 



\R{E)[' 



and 



7^ 



L^R 



mm' 



(8) 



where TZr^l represents the probability for a reflected left-handed flux from an in- 
cident hight-handed one and TZl^r that for a reflected right-handed flux from an 
incident left-handed one. From the CPT theorem we have that TZi^^ = TZr-^l and 
'R-R^l = TZl^r, where R and L label left and right antifermions, respectively [H, |^. 
We obtain for the reflection coefficients 



RiE) 
R{E) 



—e 






^11 



2iEzo ""12 



E-P 

mo 

E~P 

mo 



t^l2 



U! 



E-P 



11 



t^l2 



(9) 



with p 



E^ — rriQ. A more detailed derivation of (|^) may be found in refs. 



and [jl4|^- Now the task is to evaluate Q{z, zq) and Q{z, Zq). We consider at first the 



real mass case whence from (^) that Q{t) = Q{t). Thus, Q{z,zq) = Q{z,zq) and 
consequently TZr^l = T^l^r- Therefore, it is obvious that if the mass is real in (|^) 
no CP-violating effect exists. 

In this case Q{t) can be expressed as: Q{t) = a^E — ia2 rrio /(r), with 
m^r) = mof^r) and a are the Pauli matrices. We require f{oo) = 1 and /(r) = 0[1]. 
In what follows the characteristic energy and length, mo and ^, are used to get 
dimensionless quantities, so that eq. ([T0| ) must be written as: 



Q*{r)=a,E*-za2f{T) , (10) 

where Q*(t) = ^^ and E* = —. We do not mark the dimensionless quanti- 

^' \ / mo mo ^ 

ties with an asterisk in what follows. 

As shown in appendix A, taking into account the anti-commutation properties 
of the Pauli matrices and the definition of Q{z, zq) as a path ordered product, the 
latter can be expanded as follows 

fi(^,^o)=(Er„(E)(l)")e-^;'^^^^^) (11) 

where 

piz) = +iE'-[fiT)]Y' , (12) 

and the first terms of the expansion can be written as 

ro = i , 



T,iE) = a2 fdr p{t) /(t) e'^^^^^ />«p(«) , 

Jzo 

r,{E) = fdr rdip{T) piO f{r) /(O e-^^-/c^'^^^(^)+ 

J Zn J Zn 



^3 / „,_ „/_^/x/_^^2 



+^^ drp(r)(/(r)) 

Z J ZQ 



0-2 



dr fdi fdxp{r) piO Pix) fir) /(O fix)e^'^'''^^^'^(^)e-'"''^''"'^'^ 

Zq J Zq J Zq 



f{r){f{0?e-'-^^:'^^^^^ - /(0(/(r))'e-^^'^^/e^'^'^^(^) 



p(^)p(0 



/ J Zq J Zq 

+1 fdr p{t) {f{T)f e-2-3/>xP(x)} . (13) 

Z JZQ 

In principle, the expansion (pJ]) must be convergent as a consequence of the definition 
of VL{z, zq). Nevertheless, an explicit proof is not easy. In any case, it can be directly 



shown I^T], although tediously, that ([TTl) is a good asymptotic expansion by using a 
diagrammatic notation and certain prescriptions to obtain the mentioned expansion in 
a general way. This asymptotic character is enough for our high energy considerations. 
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In next sections we will study separately non-analytical profile functions for 
a wall with finite real thickness and analytical profile functions for a wall extending 
formally to infinity, although with a finite effective thickness because the mass goes 
exponentially to zero when r ^ — oo and to mo when r ^ +00. 



3 The non-analytical case 

We consider that the profile functions describes a finite domain wall which extends 
from r = —6^ to t = 5yj (obviously 5^ is the half of the total wall thickness). In 
general, a certain function in the region {—5yj, 5^) is matched onto /(r) = 0, 1 for 
r < —5w and r > 5^, respectively. There is no analytical function able to describe 
the wall structure from zq < —5yj to z > 5yj because at least a derivative of a certain 
order does not exist at the points r = —5^, 5^,- All the derivatives defined to the left 
of r = —5w or to the right of r = 5^ are zero, 

lim /(")(r) = lim /("^(r) = , V n > ; (14) 

nevertheless, there must be certain feg, ki that verify 

lim /('=«)(r)^0 , lim_/('=^)(r) ^0 . (15) 

The /co-derivative and the fci-derivative are not continuous functions on r, they are 
not well-defined for r = —5w,5u)i respectively and /(r) is not analytical on the real 
axis. 

In order to study this case it is convenient to introduce the new integration 
variable x = 4- and the parameter a = 2E6w. We consider the evolution inside the 
domain wall (i.e. from zq = —5^ to z = 5^). In principle we take a and E as two 
different parameters in the expansion. Following these prescriptions and expanding 
p(r) in powers of (1/-E) we can derive from ([T3| ) 



3 icT'^ax 



T,{E) = ^2e~-^^{y'^ dx F{x) e^--^^ - (i) ^ £ dx {F{x)) 

+^crs (-^)' ^ £ rfx£ dy{Fiy)rFix) e^^-} 
T^iE) = — f dx r dy F{x) F{y) e-*"^"^ e*"^'^^ + ^a / dx{F{x)) 
T^iE) = 02 e-*"^"{— j dx r dy T dt F{x) F{y) F{t) e^'^3a(a;-y+i) 



-^ , ^o. , ^y . v-A- vy;^^ "'"'"'^ - P^-^^P^-^^^ -'""'"y 

J-l J-l L 

+^ / dx{F{x)f e^"^"^} . (16) 



For convenience we define F{x) = f{6wx). By assuming that E"^ ^ a^, the exponen- 
tials containing integrals of F{x) have been expanded in ( J? )• In (|16]) the expansion 

is considered up to terms O ({1/E) ). In general, the comparison between (0) and 
(PI) gives expressions for un and u;i2 depending on the functions r„(£'), the former 
is written as an expansion for n even and the latter for n odd. Follows from ([T6|) and 
after some calculations up to order O {l/E) 



UJll 



1 /1\^ 
^^A\e) (^2,2(a)+2i?2,i(a)) + o 



E 



J,<i> 



(17) 



co'12 = —le 



^^^)-:»n(i 



(^3,3 



a) 



iRl^{a)R2,i{a)) + o 



E 



-i<f> 



where 



i?i,i(a) = a dx F{x) e* 



Ja{y-x) 



-^2,1 (a) = a dx{F{x)) 

i?2,2(a) = a^ /" dx r dy F{x) F{y) e 

RsM = «'/_/a;|_'jrfi/|_' rft F{x) F{y) F{t) e 



ia{t—y+x) 



(18) 



(f) is defined as the classical path, = J^ dr p{t), and for convenience the expansion 
of the exponentials of will be considered only at the end. 

From uu and Uu we can obtain, eq. (^, a general expansion for the reflection 
coefficient R{E). It can be seen that the numerator of (^ only contains odd powers 
of -|, while the even ones only are in the denominator. As a consequence of this fact, 
by expanding the denominator, R{E) can be written as a series of only odd powers 



in -^. It is convenient to write 



R{E) 



M<t>- 



■^^Ro{E) 



where 



-Ro(-E') = ^Cn 



n=\ 



(19) 



(20) 



It is interesting to stress that the phase factor e^*^'^"") which multiplies Rq{E) depends 
on the difference between the classical path and the classical free path neglecting the 
mass function. This phase will be of the order 0[(-^ 



1 \2l 



From eq.(|T8]), the first coefficients for (pO|) can be identified. We have 
1 



Ci(a) 



2'^ 



-iRiAa\ 



C.ia) 






(21) 



where the functions i?i,i(a), i?2,2(a) and i?2,i(a) are defined in (|T^). In order to 
express the final result in a more compact way, i?33 can be written in terms of R'^^, 
plus another term which will be canceled during the calculation. R!^^ defined as 
follows 



KM 



1 .1 

dx / dy 

-1 Jx 



dt F{x) F{y) F{t) e 



ia{t—y+x) 



(22) 



The validity of (pT]) is general provided that the profile function is strictly equal to 
for T < —5yj and 1 for r > 5^,. We will now check this result by considering a simple 
example as the step profile function, /(r) = 9{t). Then F{x) = 6{x) from x = —1 
to X = 1. Obviously the parameter 6^ lacks any physical meaning and we hope that 
the final result is independent on a. In fact, the exact refiection coefficient has been 
directly calculated for this case (see for example 0) and we have 



R{E) 



p + E 



(23) 



where p has been defined in (0). By expanding (^) in powers of -^, we obtain 



R{E) = j 

On the other hand, from ([T8|) , 
9{x) we obtain 



1/1^3 
8 V^. 



— + - — +0 



^x3- 

E 



(24) 



l|) and ( p2|) applied with the step function F{x) 



Ci{a) 



Therefore, taking into account that 
from (|T9|) and (gg) 



1/2 , 

^(l + 2za) . 



a)^ 



(25) 
in this case, it follows 



R{E) 



-^Ki) "+" 



iQ)-^a/<--'- 



^x3- 

E 



-(- 

2 V^ 



1 /I 



8 V^ 



+ 



^x3- 



(26) 



Thus, we have a positive check for (plf) . 

Now, we study the problem in a general way. By considering the arguments 
exposed in the beginning of this section, we introduce a certain function F{x) between 
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X = —1 and X = 1, in such a way that the function and all its fc-derivatives up to 
k = h—1 are matched onto F{x) = 0, 1 for x < —1 and x > 1, respectively. Therefore, 
F{x) must verify 



fW(1) =F('=)(-1) = , 0<k<h ; 

F(-1)=0 . 
With these requirements it follows from ( ]I8| ) and (^) that 



where^ 



oo -k 



Ciia) = E - A' e 



iV^k 



(27) 



(2^ 



k=h 



A, 



V'fc = TT + arctan 



(f(^)(1))' + (f('=)(-1))' - 2 cos 2a F('=)(1) fW(-1) 
i7^W(l)+i7^(fe)(_l 



1/2 



i7^(fc)(l)_F(fc)(_l 



-tan a 



To obtain (EBf), the following result has been used 



(29) 



dr F(T)e'' 



la 



E 

fc=0 



i?W( 



X) 



(30) 



Last eq. ( PPD can be usually found for any function F(x) which can be written as a 
finite power seriesp2|. Nevertheless, it is possible to prove it to be valid more generally 
for any power series provided that the r.h.s. of eq. (|30D is absolutely convergent |]2l|] . 
Moreover, by successively integrating by parts, each one of the different orders in 1/a 
appears and the asymptotic character of such series, which is the only requirement 
we need, can be proven provided that F{x) is not a quasi-periodic function (i.e. 
it cannot be written as the product of a positive function times a periodic one). 
We will therefore concentrate on non quasi-periodic profiles, although an interesting 
physical consequence of the quasi-periodic ones, which can be well understood in our 
formalism, must be shown before. The decreasing with the energy of the reflection 
coefficient and the resulting CP asymmetry (as will be seen in section 5) is due to 
the position dependent optical phase, e*"^', in our formulas. A larger energy produces 
oscillations of a minor period for the integrand in the left hand-side of eq. (|30D, in 
such a way that the integral goes to zero, decreasing as energy increases. Nevertheless, 
if the profile function is quasi-periodic positive interferences can take place. In fact, 
when the energy is taken so that both periods, the one for the oscillations due to the 
optical phase and the one of the quasi-periodic profile, are the same, the integrand 
is defined positive and a peak is expected around this value. Moreover, by taking 
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into account that a = 2E5w, this peak shifts towards larger energies when either a 
thinner wall is considered or the period of the profile decreases. For the latter case, a 
minor period will also generate a concentration of the asymmetry around the peak, 
since the effect of the oscillations due to the profile outside the energy range of the 
peak will be the same as those of the high energy optical phase. These effects were 



pointed out in ref. [^ for a particular quasi-periodic profile wall. 

Until now, E and a have been considered as different parameters in the ex- 
pansion. In fact, we expand in -^ for large E and consider a fixed. In this case 
(^) gives the leading term in the asymptotic behaviour of the reflection coefficient. 
The convergence of the sum depends in general on the convergence of the right-hand 
side of (^), and the good asymptotic behaviour of the latter guarantees that of the 
former. Nevertheless, a fixed for E large implies 5^ small and in this way the profile 
function considered is representing a thin wall, thinner for larger E. By comparing 
the reflection coefficient obtained from (^) and ( P5| ) with ([2B|), we can see that a 



power law for the fall with the energy is obtained in both cases, and in the limit 
a — !► the result for the step case is recovered. It is easy to show from ([T8|) and (p|) 
that 

C,{a) = ]^ + 0[a] , 

C,{a) = \ + 0[a] . (31) 

It results that the validity of the step profile functions demands, beyond the obvious 
condition (5^ ^ 1, the additional one a <^ 1. With the appropriate dimensions, these 
conditions are 5w <^ h/niQC and E -C Tic/5w, respectively. The latter giving the range 
of energy where the fermions are not sensitive to the real functional form of the wall. 
A more interesting situation is when 5^ is kept fixed and a increases with E. 
The expansion ( pOf ) must be rewritten in terms of E and 5yj. We obtain 

/I X fc + l 

Ro{E) = Y.'^,{5^,E)(^-^ , (32) 



where 



i^ 



V^{S^,E) = j^^A,e'^- . (33) 

The remaining dependence on E of the coefficients T>k in (|3^) is due to the factors 
g±ia -^j-^ich appear in the integrals of (p!8| ) and they are not expandable in powers of 
-^. It is worth to stress that the exponential factor e^'^*^"") in ([19|) can be expanded 
in powers of (l/E) because, in this term, the phase goes to as -E ^ cxo. This energy 
dependence of Vk does not affect the good asymptotic properties of the expansion 

(^), i.e. Vfi [-^] gives the high energy behaviour of the reflection coefficient, since 
the factors e^*°, and therefore the coefficients V^, remain bounded when E ^> oo. 
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These factors e^*" are crucial for the study of the CP-violating process, as will be 
seen in section 5, because they provide a non-trivial CP-even phase. Eqs. (|3^) and 
(P^l) are derived from (PD|) and (P^D provided that the expansion in eq. (^) keeps its 
asymptotic character as the dependence on E of the coefficients Cn{a) through the 
parameter a is considered. The latter happens when 

imi ^ , = (34) 

for all odd k. The observance of this condition can be proven although it is not an 
immediate consequence of the asymptotic character of the expansion in (|ll|) because 
the profile function /(r) is not zero outside the domain wall : /(r) = 1 for t > 6w 
It follows from (pTf) that (|34D is satisfied for k = 1. A general proof may be found in 
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Therefore, for the functions introduced in this section as mathematical rep- 
resentation of a strictly finite physical wall, we obtain a fall with the energy as the 
power —{h+ 1), where h is the order of the first derivative which is not continu- 
ous (or does not exist) at any of the two point r = 0, 1, or at both. This result 
can be easily generalized in order to assert that the dominant high energy behaviour, 
when the profile function is h times derivable on the real axis, its h-derivative, being 
discontinuous, is ~ 1/E^^^ as in eq. (|5^. This is the main result of this section. 



4 Analytic profile functions 



We consider now the case where the profile function is analytic on the real axis. By 
considering the arguments we presented in the previous section it follows immediately 
that an analytic function cannot describe a wall profile which is constant outside the 
domain wall. Therefore, these profile functions can only be asymptotically equal to 
0, 1, as r ^ — oo, +oo, respectively. Thus, the wall extends formally from r = — cxd 
to r = -|-oo. Nevertheless, this wall will be characterized in general by an effective 
thickness, a, which allows to define a certain domain wall. The particular criterion 
considered in order to define a is not important. For instance, |/(t) — 1| < 0.1 for 
T > a and |/(t)| < 0.1 for r < —a. However, in general we can write /(r) = F f M, 
where we will assume 

F{x) ~ e^^- , F{x) ~ 1 - e-^^+ ; (35) 

being C+ > 0, C_ > 0, parameters depending on the particular profile wall. The 
asymptotic functional behaviour assumed in eq. (^) is in principle the simplest that 
can be demanded to a profile function verifying the initial asymptotic conditions as 
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r — > — cx), +00 introduced above. Furthermore, it is not easy to find functional 
beliaviour otlier tlian tliat of (|35D, with good analytic properties on the real axis, 
and, on the other hand, a naive derivation of the wall profile from the corresponding 
equation of motion as done in [Q agrees with this assumption^. We consider the 
evolution from r = —6^ to r = 6^, but this parameter is now defined in such a way 
that f{5wx) = F{Xx), where A = ^ verifies 1 ^ A ^ ■^. Thus, analogously to the 
previous section, it follows from (|TTD and (|l^) 



f 1 a /■! 

Q{Xa, -Xa) = < 1 + -^^0-2- / dx F{Xx) e~ 
E 2 J -1 



■io"3a(l— x) 



+ 



E 



^10-30 



(36) 



where all the previous definitions are kept. 

Again, we define R{E) = e^'^^'"-^ Rq{E) and obtain 



UE) = \ (i 



la 



dx F{Xx)e' 



+ 



E 



(37) 



In order to derive ( P7[ ) from (PB]), the first of the equations (|^) must be used, and 
therefore the profile function is modified outside the region (—6^,6^) to describe 
a situation with zero and constant mass for r < —6^, t > 6^. This technical 
assumption, necessary to identify the incident, reflected and transmitted particle 
fluxes outside the region considered, introduces non-negligible effects for the high 
energy behaviour which must be separated from the real ones caused by the evolution 
through the wall region considered. In fact, the numerical calculation performed in 
requires this technical assumption and as a consequence the authors consider a 
wall strictly located in the region between z = and z = zq despite they work with 
an analytical ansatz which goes exponentially to zero and one when z —>■ —00, +00, 
respectively. However, this modification has no practical consequence except for the 
high energy behaviour. It should be noted that the hypothesis E ^ 6^ is necessary to 
obtain (^), therefore we cannot take the limit 5^; — >■ 00 directly in order to calculate 
the reflection coefficient. When the profile function varies from r = —00 to r = +00, 
the coefficients must be formally defined by using the asymptotic fiux at r = +00, 
—00. Obviously, it should be done before obtaining eq. (|36D , by considering the 
evolution of the fermions from r = —00 to r = +00 , without the introduction of 
the parameter 6^ and without any additional condition for the energy beyond the 
initial one, E ^ 1. Nevertheless, in this case we cannot expand the exponentials 



of ([T3| ) to derive an analogous expression to (|T6D, the integrals we obtained for the 
expansion of the reflection coefficient are then much more complex and practically 
it is not possible to analyze them in a general way. This is why we approach the 
problem in a different way: we study in practice the evolution from —S^] to Sw in 
order to find the total contribution from —00 to +cx3 and other terms depending on 
6u], which must be in principle negligible by considering the requirement ( ^5]) for the 
profile function. 
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By rewriting appropriately the integral in ( P7| ) as the sum of two terms, the 
first giving the contribution of the whole integration domain and the second due to 
contributions of a; ^ [—1,1], we obtain after some calculation, for the leading term. 



Ri^\E) = -I (-) I {^m /_ J dx F{Xx) e^^ )^} + (-) R^.c{a) (St 



where, from (|35D and A ^ 1 



RN-c{a) = +2 \ 



-AC+ ia ^ I „-AC_ „-ia '^ 



e ^ e 



ta — AG+ la + AG_ I 



It can be seen that the first term of the right-hand side of (^) does not depend 
on Su, (Although it seems to have a dependence on 5^ through a, if we change the 
integration variable x by a;' = Xx, we see that this dependence disappears) . Now, we 
analyze its functional dependence on the energy and then we turn to the question 
whether Rn^c{(i) is negligible with regard to this term. 

If the Laurent expansion for F{z)^ the analytic complex extension of F{x), in 

+ 00 

the pole of order z/^, z = Zj, is J2 (^ni^ ~ ^j)"' the following result can be shown by 

n=—i' 

applying the Cauchy theorem 

lim / dx F{Xx) e(*»-^)" = 27rz Y e"^^^ e^^^' Y aL^-^ . (40) 

e^o+J-^ ^ ' ^tl ti A"(n-1)! 

Where zi = xi+ iyi, Z2 = X2 + iy2, ■■■, zn = x^ + iyN, are all the poles of F{z) with 
positive imaginary part we have picked when the integration contour is adequately 
closed. As the profile function, F{x), is analytic on the real axis, it is obvious that 
yj 7^ and therefore we obtain an exponential dependence on a. Nevertheless, in this 
case the parameter a has no physical meaning because 6^ is arbitrary. The distance 
6u] is taken as several times the parameter a in order to apply eq. (^). Moreover, a 
is defined after an arbitrary criterion. The right-hand side of eq. (^O]) may however 
be rewritten as follows 

where the dependence on 6^, through the parameter a, is replaced by the one on a. 
The final result will be expressed, (47), in terms of the singularities of f{z), instead 
F{z), eliminating any dependence on our arbitrary parameter a |jl|. 

Now, the discussion of (^T]) is very similar to the previous one for the non- 
analytic case. If the range of energy and the effective thickness, a, allows to consider 
that the factor Ea is not large, then the term ^R^-c is negligible with regard to the 
first one of eq. (|38|) , 

aE e-^^^y^ > e-^^+ , e"^^" ; (42) 
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since A is large^ the latter gives then the asymptotic behaviour of the reflection 
coefficient. Nevertheless, by considering the energy dependence of (|l|) and Rn_c{(^), 
it is obvious that for large enough energy the first term of (PB| ) is smaller than ^R^-c 



since p^) is not satisfied. This result is at first surprising and seems to disagree 
with our requirements about the profile function. But, now we must consider the 
above mentioned technical assumption since as a consequence of that, the profile 
function we are studying in practice is non-continuous at r = — 5^, and r = Sw 
Therefore, a contribution of the order -| to some power follows from the results of 
the previous section which obviously must depend on 6^. It is not easy to obtain 
exactly these contributions because the result (^) cannot be in general applied to 
these profile functions, except in the large a limit where the asymptotic behaviour 
can be determined. It can be proven that the asymptotic behaviour for large E, and 
therefore a, of the integral in the left-hand side of eq. ( pOD is indeed given by the first 
term in the right-hand side of eq. (PU]), in a completely general way for any function 
F{x). Thus, the contribution of the non-continuities can be identified. Moreover, if 
we take the limit a — ;► then Rn-c ~~^ 0, which is the behaviour we hope for the 
contribution of non-continuities: the profile function gives the step function, in this 
limit, and it can also be proven that the contribution of the right-hand side's first 
term of eq. ( pSD is exactly the refiection coefficient for the step case |]2^ . The energy 
dependence of the term ^R^-c is explained in this way, as well its surprising large 
energy dominance. However we are interested in the high energy contribution of the 
total evolution from r = — cxo to r = -(-oo, and all the previous discussions allow to 



assume that it is given by the first term of eq. (|38|) . By using eq. (|4l|) we obtain for 
this high energy behaviour of the refiection coefficient 

RiE) = 2na ^ e'^^^^^^ e^^^-^^^ ^ ai /f _; . (43) 

j=l n=l ('^ -'^)- 

Besides, if we are in the range of energy where Ea is large, the leading contri- 
bution is given by the pole of F{z) with the smallest imaginary part, i.e. the closest 
to the real axis. Calling z^ = Xk + iyk this prevailing pole we obtain 

RiE) = 27rae~2^'^^^ a^^ ^J^l^i^ ''''''''''' ' ^^^^ 

In this case there is no proof of a result analogous to (0) (even the particular checks 
are difficult). However, by studying directly ( [TB| ) as well as the prescription to obtain 



the coefficients of the general expansion (O) which is given in ref. [21], it can be 



reasonably guessed that the asymptotic character of eq. (PD|) is given by 



hm ^^ — r <a^ . (45) 
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Looking at ( [TB| ) and taking into account the expansion of eq. (|T^ for p(r) in powers 
of (1/-E), it is easy to see that T^i^E) contains terms of the order E'^a^, because the 
first of the three integrals of T^lE) in eq. (|T^) appears multiphed by {crE)^ when 
the change of integration variable, a; = -, is performed after the introduction of the 
function F{t) = /(err). If we conjecture that this result can be generalized to Tn{E) 
containing terms of the order E'^a"' (this conjecture is supported by the study of the 
general expansion given by eq. ([Ill) in ref. [^), we get the ratio given by eq. (|45|) . 



If the particular profile function considered generates an expansion for the reflection 
coefficient which verifies (^Sj), then (^31) and (^4]) will be valid in any range of energy 
satisfying E ^ 1 and o" ^ 1 (obviously (0) is only valid for aE large). As can be 
seen above, these conditions are suitable for studying quarks propagating through a 
wall except for the case of the top. 

We worked by convenience with the function F{z), obtained through the 
change of variables performed at the beginning of the present section. Let us now 
return to the physical variables and express the final result as a function of the sin- 
gularities of f{z) instead of F{z) in order to avoid the apparent arbitrariness arising 

from the dependence on a. If we consider the Laurent expansion J2 (^ni^ ~ -^j)"? 

n=-Uj 

for f{z) in the pole of the order Uj, z = Zj and take into account that f{z) = E{^), 
we obtain 



Zj = azj , 

a^ = ^ • (46) 

a 



Eq. (M) can thus be re-written as follows 



R{E) = 27ie-^^y'^ a% i^l3ll g^^^^ ; (47) 

[J^k - 1)! 

where the former condition aE 3> 1 must be re-formulated by requiring that e^^^'^^'yu) 
may be neglected, y^ (y^) being the imaginary parts of the (next to) closest pole to 
the real axis. It is also interesting to notice that the wall thickness we introduced 
through the parameter a is now incorporated without arbitrariness into the profile 
function. From eq. (^61), the connection between this thickness and the distance 
between the poles of f{z) is apparent. 

Now, we check the present results by using the particular analytic solution 
obtained in references [^ and [|T^ for the ansatz 

/W = i±^4<^ . (48) 

For this profile function it follows from (^) that 

R(E) = 7r .. X , (49) 

^ ' 2sinh(7rEa) ' ^ ' 
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which agrees with the resuh of [|T^ and [|I^ in the above mentioned conditions. 



Therefore, with this positive check, we conclude that (^31) in general and (|J) 
for Ea large give the high energy behaviour of the reflection coefficient, for the 
analytic profile functions: at high energy the reflection coefficient decreases exponen- 
tially, the coefficient in the exponential being the smallest positive imaginary part of 
the complexified profile function. This is the main result of this section. 

It is interesting to emphasize the exponential dependence on the energy in 
this case instead the power dependence in the previous one. Apparently, there is no 
connection between power law energy dependence derived in in the preceding section 
for non-analytic profiles, and the exponential law for analytic profiles derived in the 
present section. Nevertheless it can be proven that the results for non-analytic profiles 
can be understood through the properties in the complex plane of analytic extensions 
for these profile functions by using eq. 



5 CP-violating wall profiles 



As pointed out in the introduction, our effort to study CP-violating effects concen- 
trates on CP violation at the tree level due to the expectation value of the Higgs fields 
beyond the Standard Model. The time independent Dirac Hamiltonian incorporating 
CP-violation can be written as follows 

H = dp + (3mR{z) + iP'j^miiz) . (50) 

As stated in ref. |0, m^^z) and mj{z) are real because of the hermiticity of the 
Hamiltonian. The CP violation comes from the term oc mj{z), when it cannot be 
rotated away through a chiral rotation. Therefore, we must go beyond the SM in 
order to find CP violation at tree level. In the SM, with one Higgs field mj{z)/mFt{z) 
is a constant independent of z and therefore this source of CP violation can be rotated 
away. In the other models, with more than one Higgs field, mj{z)/mfi{z) depends in 
general on z and the matrix mj{z) will remain. 

Following the notations and conventions of section 2. the time independent 
Dirac equation can be written as in eqs. (^ and (Q), where m{z) = mii{z) + imj{z). 
The CP-violating effects are so incorporated by assigning a complex mass to the 
fermion inside the domain wall. 

Analogously to (p!0D , we define mji{T) = rnof^r) and mj{z) = emQgij), where 
/(r) and g{T) are real functions of order 1 which characterize the CP-violating profile, 
and e is a parameter associated to CP-violation. We express 

Q(r) = a^E - i {cr2f{T) - eaig{T)) , 
Q{T)=asE-^{a2f{r)+ea^g{r)) . (51) 
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We keep the same diinensionless quantities, the same definitions of the previous 
sections and assume that g^r) goes exponentially to zero at both, r — >• +00, r -^ —00, 
for an analytic profile or that it it is strictly zero outside the domain wall for a non- 
analytic one. This choice is always possible by an appropriate chiral rotation which 
rotates away the imaginary part r?7,/(+oo) , i.e. which corresponds to a real mass in 
the broken phase. We have checked that the same result is obtained with another 
choice. 

Thus, the high energy leading term for the refiection coefficient can be ex- 
pressed in a general way as 

where G{x) is defined as F{x) in the previous section, x being the same dimensionless 
variable which was then introduced, and A = 1, 6yj/a for non-analytic, analytic 
profiles, respectively. The sign — is for R{E) and + for R{E). 

For both, analytic and non-analytic profiles, the difference TZr-^l-T^l^r, which 
is a quantitative measure of the CP-violating effects induced by the wall, can be 
computed. For the sake of simplicity we concentrate on the non-analytic ones: 

7JH..-7J..H = -4.(i)''"'"ji^sm(A^) , (53) 

where the ^-derivative of G{x) has been considered as the first non-continuous one, 
where Bi is defined from G{x) as A^ in ( pQ] ) from F{x), and where 



Aip = ae-^h + {i-h)^ , (54) 

with (fh defined in ([29| ) from F{x) and ag analogously from G{x). When gi^r) vanishes, 
Bi = and as expected the r.h.s of ( [^ ) vanishes. 

The former requirement, g{T) vanishing outside the domain wall, is imposed 
as above mentioned in order to work in the physical basis, in which the masses are 
real. For technical reasons, we will also consider the chirally rotated convention in 
which g{T) = 1 in the broken phase outside the domain wall, i.e. mi = erriQ. 

One can easily check from (^Tj) that the imaginary mass in the broken phase 
will be removed by the chiral rotation exp(ir775) with e = tan(2?7), and the real 
mass in this physical region will become m = mo|vT+T^|. As a matter of the 
fact, in that region the eigenvalues of the matrix Q{z) and Q{z) will be ±p, where 
p = IVE"^ — m?\^ m obviously being the physical mass. Therefore, the dimensionless 
quantities we introduced in previous sections should be now built by using m instead 
niQ. Then, eq. (^) must be replaced by 

(-a!!''! I - - f4) {^-^^ -^^Cdx FiXx) e-^ (55) 

\R^'{E)j 2\EJ\ ^/TT^ ./T^^J^i ^ ^ ^ ^ 
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ae 



vm^ 



1 

dx G{\x) e*' 



and it follows that 



nn^L-n,^n = -Y^,[-^) j^^-y^,MA^) , (56) 



instead of eq. (|53|) . By comparing eqs. (^31) and (|56|) it is very easy to see that all 
physical consequences and in particular the discussion which follows will be valid for 
both cases. For example, the vanishing of the r.h.s. of (|56D when no CP violation 



is present is now due to Aip = 0. An amusing remark is that eq. (|56D shows that 
the highest chiral asymmetry originated by the CP- violating wall profile appears for 
the case e = 1. In other words, by taking into account that as well /(r) as g^r) are 
required to be 0[1], the most effective situation in what concerns the production of a 
net baryon number left in the broken phase for the high energy range arises then from 
requiring that real and pure imaginary mass terms are of the same order. In fact, Eq. 
(^) remains unchanged if we make £—>•-, this symmetric behaviour for £ > 1 and 
e < 1 being coherent with the fact that real and pure imaginary masses in the broken 
phase can be turned into each other through a chiral rotation. As already stated, 
the imaginary mass in the broken phase can be rotated away by a chiral rotation 
to return from the latter convention, g{oo) = 1, to the former one in which g{T) 
vanishes outside the domain wall. The chiral rotation implies a redefinition of the 
new real and pure imaginary functions, J^(t) and G{t), as a function of the old ones, 
/(r) and 5'(t). If we then apply the new functions, J^{t) and Q{t), which behave 
now appropriately outside the domain wall, to eq. (|5^ ) and obtain the first non-null 
contribution to the chiral asymmetry, eq. (|56D will be re-obtained. This expected 
result is not immediate to obtain explicitly, because the contribution to the chiral 
symmetry given by eq.(p3D, which is in general non-null, gives zero for J-'{t) and Q{t) 
ii h ^ i. Therefore upper order terms must be considered to re-obtain eq. (|5^. 

As well known, a CP- violating effect is always generated by the interference 
between a CP-even phase and a CP-odd one. In fact, it can be seen from eqs. (E 



and (p3) that a non-zero leading term first needs the CP-odd phase arising from 



the reversing of the sign for the pure imaginary mass term in (52), but also a non-zero 
angle Ay?, i.e. a net CP-even phase. Notice that working at high energy {E ^ 1) we 
are far off the total reflection domain {E < 1), and that due to the thickness of the 
wall nothing prevents a CP-even phase. In fact the CP-even phase is generated by the 
interference of the phase of the mass term mj{z)/mfi{z) at different values of z. In 
the step approximation, it is always possible to make the mass term real by rotating 
away the phase through a chiral rotation, but if we have a z-dependent complex 
mass term inside a certain region, this rotation is not possible in general and the 
imaginary term generally remains. Tree level effects can thus appear. Nevertheless, 
if the z-dependence of the imaginary term is proportional to the real one, we can, 
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again through a chiral rotation, turn the complex mass into reaL When we take 
the convention g{+oo) = this situation corresponds to g{z) = for all z. This is 
stressed in references ||^ and ||T^ and it is obvious from (PD| ) that in this case h = i, 
iph = ae and therefore Aip = 0, which lead to TZr^l-T^l^r = 0. Furthermore, the 
leading term will be zero if we assume a weaker condition than the proportionality 
of real and pure imaginary terms. 

Assuming F{x) = C + Fo{x) and G{x) = K + Go{x), where C, K are simple 
constants and Fq{—x) = —Fq{x), Go{—x) = —Go{—x), we obtain 

fW(1) = (-l)'*+i f(^)(-i) , 

GW(1) = (-1)^+1 F(^)(-1) ; (57) 



ioT h, i > 0. It follows from (|57D and (|5) that under these conditions 



Aip = n7i , (58) 

where n is an integer. Analogously we obtain (^) if both functions have even sym- 
metry with regard to the axis x = (It is obvious that our requirement F{+oo) = 1 
forbids an even symmetry for F{x), but our results can be applied to general situa- 
tions where F[+oo) = 0). We conclude therefore that this weaker condition implies 
a null CP- violating leading term as it follows from (0). This does not lead to a 
strict suppression of CP-violating effects as the stronger requirement of proportion- 
ality does, but it can be stated that CP-violation will be of a lower order in the high 
energy range. Whether these symmetry conditions involve any similar consequence 
in a lower energy range is an interesting question, but it is out of the scope of the 
present work. 

It is also interesting to stress that the most favourable situation is given by 
the opposite symmetry properties of F{x) and G{x) with regard to the axis a; = 0. 
If we assume that the functions Fq{x) and Go{x), above defined, are even and odd 
symmetric, respectively, or vice versa, we immediately obtain 

Av, = ^vr , (59) 

n being an integer. The highest value for sin(A(y9) is thus given by imposing these 
symmetry requirements. 

An expression analogous to ( ^31) can be found for the analytic case by using 
(0), and a similar discussion may be done. The crucial point for such a discussion is 
that the integral /^ dxF{Xx)e^^"'~'^^^, which is obtained from f_^ dxF{\x)e^"'^, is pure 
imaginary if F{x) differs by any constant from an odd symmetric function, and real if 
it is even symmetric. It is very easy to see that, for instance, the analytic ansatz (^) 
is odd symmetric with regard to the axis x = up to a constant, its contribution to 
the reflection coefficient has indeed no CP even phase factor (eq. (p9D). If a function 
F{x) for the real mass term and a function G{x) for the imaginary one have both 
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odd symmetry, the angle Ayj vanishes. No CP-even phase is obtained and no CP- 
violating effects at tree-level and at first order in the high energy limit arise from this 
case. Nevertheless, if a function G{x) defined as —^, which is even symmetric, is 
taken, a net CP-even phase and therefore CP-violating effects at the first order are 
obtained. This simple ansatz for G{x) illustrates our statement about the symmetry 
properties. 

In ref. [^ numerical calculations are performed by considering the CP- 



violating term as a perturbation to the CP conserved Dirac equation, for the partic- 
ular real profile function above mentioned, F{x) = 1/2(1 + tanhx), and for the two 
following pure imaginary ones 

G2{x) = {F{x)f . (60) 

The quantity A'"^ = IZr^l — TZl^r for both particular cases can be immediately 
obtained in the high energy range from the general results given by eqs. (^31) and 



A^^ = -47r2£a3Ee-2-^'^ . (61) 

These results are valid obviously for E ^ 1 and a ^ 1 as we explained in section 4, 
but the requirement Ea ^ 1 is also necessary to isolate the contribution of the leading 
singularity EM. Nevertheless, the contribution for all the poles can be resummed and 



we obtain, for instance, in the first of the two former cases 



Ar = 2.a^E . ; , (62) 

which is a valid result for any Ecr. 

The energy range we are interested in is not explored in the numerical results 
given in ref. [0, although a few results for E = 5.0 are presented amusingly in 
relatively good agreement with the ones we obtain from our analytic expression^. 
We also checked in section 4 that the unperturbed refiection coefficient which is 
obtained in that work for the above expressed particular ansatz completely agrees 
in the appropriate limit with our general result for a non CP-violating wall profile. 
Furthermore, we obtain the ratio Af^/A^^ = —2 which is observed for the results 



given by Funakubo et al. |]19|, even for E ^ 1. 

Several interesting remarks must be finally done. The function G2 (a^) = [-^(a;)]^ 
does not present (up to a constant) any well defined symmetry with respect to the axis 
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X = 0, while Gi[x] = —^ is even symmetric; a higher CP asymmetry for the latter 
than for the former is thus not surprising in view of our conclusions about symmetry. 
Moreover, the dependence on a of eq. (p^ also agrees with the well-known fact that 
no observable CP-asymmetry can appear in the thin wall case at tree level, i.e. in 
the limit a ^ 0. Funakubo et al state that the degree of decreasing of A*^^ seems 
to be larger as the wall thickness increases after studying two particular cases for 
several values of the wall thickness. Regarding eq. (^) we can state a more precise 
conclusion for any non-analytic wall profile in the high energy range which confirms 
the trend of the CP asymmetry to decrease when thickness increases pointed out for 
those authors. The quantity A'-^'^ generally is an oscillating function on the thickness 
parameter because the angle Ay) depends on it, but those oscillations are modulated 
by a factor decreasing with larger thickness as eq. ( [53| ) shows (A*^^ as a function of 
the thickness has in general an infinite number of zeros, the first of them being in 
fact as the thickness is zero). The functional behaviour of A"-^"^ given in eqs. (pH) 
and (|62D for analytic wall profiles seems also confirm that trend, but in that case our 
results are only valid for a <^ 1. 



6 Summary and conclusions 



Rigorous and completely general conclusions about the behaviour of fermions hitting 
a wall with an arbitrary profile are very difficult. Obviously, the fermion scattering 
is determined by the particular profile function which characterizes the wall we are 
treating. The crucial point in order to get general answers in this problem is to show 
relations between the functional properties of the wall profile and the refiection co- 
efficients, etc, which characterize the behaviour of the fermions hitting the wall. A 
prescription to solve the Dirac equation in the presence of CP-violating electroweak 
bubble wall is presented in [^]. Nevertheless, the results depend on certain functions 
of the energy which must be obtained solving a second-order differential equation. 
This prescription, applied to the particular ansatz studied in ]TB|], leads to a hyper- 
geometric differential equation, but in general it is a difficult problem which does 
not allow to relate in a direct and useful way the wall profile function to the fermion 
behaviour. 

A formalism to study the scattering of the high energy fermions in the pres- 
ence of a general wall is proposed in the present work. High energy fermions mean 
E/rriQ ^1. If we consider that the temperature during the supposedly electroweak 



first order phase transition is about 100 GeV |T6|, the Boltzmann thermal distribu- 
tion gives an average energy for the quarks which verifies ^ = O[10~^] except for the 
top. Thus, the range of energy considered is interesting in the cosmological problem, 
although for the top we only study a high energy tail of the thermal distribution. 
The quantum corrections to the classical behaviour expected for this range have been 
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obtained for two different types of mathematical representations of the wall. The first 
describes walls with profile functions which differ from the asymptotic values only in 
a finite domain. Whether they are realistic is a difficult question. They are modeling 
situations involving different length scales of variation for the profile function and 
its derivatives, i.e. so that to a given accuracy, some derivatives can be considered 
as discontinuous. In the second case, when we assume that the wall profile function 
differs from the asymptotic values up to infinity, edge effects related with different 
length scales of variation are neglected (quantum, thermal and statistics fluctuation 
are not considered, for example). However, we study the quantum corrections in 
both cases and show how the high energy fermions are sensitive to these facts. As 
can be expected, the quantum corrections are more important as the profile becomes 
sharper. In fact, if the h-th derivative of the profile function is the first non contin- 
uous one, i.e. it presents a certain step, the high energy behaviour of the reflection 
coefficient is characterized by the following power law 

R{E) oc moc\hc/6^)\l/E)''+^ , (63) 

against the exponential behaviour 

R{E)<xexp[-yiE/{hc)] , (64) 

given by a profile function analytic on the real axis. The physical dimensions are re- 
stored in eqs. ( |55[ ) and (0). S^ in the first equation is the wall thickness parameter 
defined in section 3 and the coefficient y'j^ in the latter one is the smallest positive 
imaginary part among the poles of the analytic complex extension of the profile func- 
tion /(r). It worth to stress that j/^ has the dimension of a length, since the variable 
r for /(r) characterizes the position in the normal axis to the planar wall. Notice also 
that, while in the former case the reflection coefficient is proportional to some power 
in h as expected for a quantum effect, in the latter case the reflection coefficient is 
not even analytic in h, similar to a tunnelling effect (for example, the term given by 
eq. (|6^) has an analogous dependence on h to that of the so-called Gamow factor, 
g-27rzz e /Uv ^ which gives the dominant behaviour for processes involving nuclei with 
charge Z'e and velocity v tunneling through the Coulomb barrier due to a charge 
Ze[^). 

Finally, CP-violating effects have been incorporated by introducing a complex 
mass, as we explained in section 5. The reflection coefficients, R{E) and R{E), are 
expressed as a function of both, real, chir ally-even, and pure imaginary, chirally odd 
mass functions, f{z) and g{z). We first perform a chiral rotation to have a real mass 
for z ^ oo, deeply in the broken phase. 

If g{z) does not identically vanish, the phase of the mass depends on the chi- 
rality of the incoming flux. The resulting chiral asymmetry leads to a CP-asymmetry 
through the CPT theorem. This CP-violating term introduced in the Dirac Hamil- 
tonian produces a CP-odd phase. This CP-odd phase is not enough to generate an 
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observable CP asymmetry. It is known that a CP-even phase should exist and in- 
terfere with the CP-odd one. That effect is clearly shown in our formalism through 
the dependence of the relevant quantity A'-^^ = TZr^l — T^l^r on the angle Aip. A 
non-zero angle is due to the CP-even phases arising from the integration of the real 
and pure imaginary terms of the profile function, although the asymptotic regime 
considered here is far off the total reflection. We found that the respective symmetry 
properties of the f{z) and g{z) play an important role, and we obtain the following 
general conclusions: i) CP-even phases and therefore CP-violating observable effects 
at tree level will be found in the thick wall case, provided that the symmetry proper- 
ties are not the same for both, f{z) and g{z). ii) CP-violating effects will be of lower 
order in 1/E (if not vanishing) when the symmetry properties are the same for both 
functions, iii) the most favourable situation in order to generate CP-asymmetry in 
this high energy range is given by a profile where real and pure imaginary functions 
present opposite symmetries. 

It is finally very important to remark that the axis z = is an arbitrary 
reference for the symmetry properties of functions. As expected our results are ob- 
viously independent of any shift on z, the physical conclusions thus depending on 
the respective symmetry properties of the real and imaginary mass terms. It is also 
important to notice that the conclusion about symmetry properties are valid even if 
the symmetry properties appear after a chiral rotation which in general leads to a 
non vanishing g{+oo). 
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A The integral expansion 



In this appendix we derive the integral expansion given in equations (|TTp and (|I^). 

The starting point is the definition given in (^) for Q{z, zq) which leads to the 
solutions of the time- independent Dirac equation as shown in @ . We had 

n{z,zo)=V e^^o''^'^^^^ , (65) 

where Q{t) is expressed by means of the usual Pauli's matrices in (|10D- By considering 
a certain path partition {zq, zi, ..., z^^i, zn, z), we can write 
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V e •'''0 



iT dr Q{t) _, i r dr Q(t) i T^ dr Q{t) i p dr Q{t) 



(66) 



If we take Q{t) as constant in each interval {zj^ ^j+i), defining for each one fj as the 
following integral average value 



1 r^j+i 
/, = ^/ drfir) 



(67) 



-'J "^3 



where Aj = Zj+i — Zj, it can be written by using ([T0|) 

Taking into account the following result for two operators verifying {A, B} = and 



nA+/3B 



cosh 



(a' + (3' 



1/2 



aA + pB . , 
H 7777 smh 



(a2 + /52)i/2 
which can be easily proven, it follows from (^) that: 



{a' + (3' 



1/2' 



(69) 






E 



sin(pjA 



j^jj 



(70) 



where pj = +{E'^ — ffY^'^. In order to obtain (^) — and — have been expanded in 

powers of -|. Eq. (|^) is also valid for any set of unitary operators anticommuting 
with each other and consequently analogous result to the following can be obtained 
for the complex mass case. 

By substituting the path ordered products of the right-hand in ( |H^ ) by (|7D|), 
multiplying and reordering the terms, we obtain 



i°'3j2Pi'^J 1 ^ -^'^s J2 Pi^3 ^'^sY.Pi^i 

VL{z,Zo)=e ^=0 +-a2^sin(pfeA,,)/fce ^=^+' e ^=o 



E 



k=0 



+'!)'< 



N fe-1 771-1 

N k-1 ia-3 Y^ pjAj -ia3 ^ pjAj ia^ ^ pjAj 

Y, Y. sin(pfcAfc) sin(pMA^)/fc U e ^=*+^ e ^=™+i e ^=o 

fc=0 m=0 



.03 

2 



AT 



iV 

+i^YMPk^k)f^e -« 

fc=0 
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3 ( N fc-1 m-1 



(-1 s d r iv K—L m— i 

p ) ^2 S XI 51 H sm(pfcAfc) sin(pmA^) sm(p„A„)/fc /„ /„ 
^"^ [ fc=0 "1=0 n=0 

N fe — 1 m— 1 n — 1 

X e J = fc + 1 g j — m + 1 g j = n + l g j=0 

N fe-1 



Af A;— 1 — iiTs y^ Pj^j ^^3 y2 Pj^i 



2 



A;=0 m=0 

JV m-1 



^^^ Af fc-l -i<j3 ^ pjAj ias J2 PAi 

— TtY^Y. sin(pfcAfc) sm{p^Am)fkfm e ^=™+i e ^=o 



2 



fc=Om=0 

AT fc-1 

X Af -io-3 J2 Po^i ^"^ Z) P^^J 

+ - X sin(p,Afc)/3 e ^-^+1 e ^-» 

^ fc=0 



+ ...,(71) 



where the sum in the exponential is taken to be zero if the lower index is bigger 

than the upper. The general result for two operators A, B, which verify {A, B} = 0, 

e^A = Ae~^ has been used. We know that sin(pjAj) = pjAj + 0[A'j] and in the limit 

Aj -^ 0, J2Aj ^ J dr with fj -^ f{r) and pj -^ p(r). These replacements can be 

j 
understood if we consider the definition of f) as an integral average value, assuming 

that fj = f{Tj) with Tj G [zj,Zj+i], provided that /(r) is continuous. Following the 

last prescriptions, Eq. (0) and (|13D can be immediately derived from ([7l|). 



B General Orthogonality properties 

In the present appendix we consider in a general way the problem of the orthogonality 
of different Dirac Hamiltonian eigenstates in the presence of a wall. We keep the same 
framework presented in section 2, assuming a z-dependent mass in order to study the 
one flavour quark propagation problem. Let {E,px,Py,Pz) be the four-moment of 
the quark. For a static wall the particle energy, E, is conserved. In the unbroken 
phase E = Pz- The system is symmetric with respect to rotations around the z-axis 
that implies conservation of total angular momentum in the z direction, J^. It is 
also invariant under Lorentz boost parallel to the x — y plane as we stressed above. 
Thus, in order to label the Hamiltonian eigenstates the conserved quantum numbers 
E,Px,Py and jz (the eigenvalue of Jz) can be used. If we boost the reference frame to 
obtain Px = Py = 0, the helicity states of the incoming plane waves in the unbroken 
phase correspond to eigenstates of Jz and this helicity may be also used to label them. 
We define the eigenstates tpn by 

H^n = {-lyEni^n, (72) 
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where we choose n = Px,Py,jz, En,r. r = 2, 1 label the positive and negative energy 
states, respectively. In [5], for the Hamiltonian 

H = ap + (3moe{z), (73) 

it is stressed that the dimension of the eigenspace is two outside the total reflection 
energy range. In this case the following two set of eigenstates can be built 



Hi + RWifne'^^^^oiz), h = j,. (74) 

where Uh{i?"''^) is a solution of the Dirac equation in the unbroken phase and the 
reflection matrix R is given by 

R = ^; (75) 

and 



^ni^) = \Fr {{usif'Y'^'''' + Jusif'Y'^'''') e{z) 

\ ; r)' L \ / 

V rz 

+ {l + J)us{f'y^^''"%-z),] (76) 

where J is the reflection matrix when the particle is coming from the broken phase, 
given by 

l + j = ^(l + i?), (77) 

Pz 

and s is a spin index dependent on jz, such that {1 + J)us{j}'^) = UhijP^^) is a massless 
spinor with helicity h = —jz', Us{f^^) satisfies the Dirac equation in the broken phase. 
The incoming, outgoing, transmitted and broken incoming four-moment are defined 
as follows 



P™"= {E,p.^,Py,Pz), 
P°"* = {E,P^,Py,-Pz), 

P^'' = {E,Px,Py,P'z), 

p'- = {E,p,,Py,-{pzr). (78) 

The wave functions ip^^'^ and ip'^ are orthogonal and allow to build an or- 
thonormal basis of the eigenspace beside the analogous antiparticle wave functions 
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(for r = 1 the spinor solution of the Dirac equation, Uh, must be replaced by the 
negative energy ones, Vh- In what follows our conclusions are valid as well for the 
r = 1 and r = 2 eigenstates). For general walls the argument about the dimension 
of the eigenspace can be generalized and the wave functions obtained to verify 



^" ^^ ^1(1 + i?)n,(j?-)e^P"'--^ z » 6^- ^""^ 



where R is the reflection matrix in each particular case and 6^ defines the character- 
istic wall thickness as can be seen above. Analogously for ip^. Thus, iplJ^'^ describes in 
the general case an incident plane wave far from the wall coming from the unbroken 
phase, bouncing on the wall, and generating a reflected plane wave in the unbroken 
phase and a transmitted one in the broken phase, ijj'^ describes the same process 
although coming from the broken phase. In [1T3, wave functions like i/j^^'^ and ijj'^ 



are obtained for the particular wall presented in section 5. In this work the authors 
calculate explicitly the overlap integral of these wave functions obtaining a non-zero 
result, concluding that they are not orthogonal. Nevertheless, a general argument 
allows to assert that ^^""^ and ip'^ must be orthogonal as we will see. In |^ the error 



in the overlap integral calculation of ref. [18| is shown and it is explicitly proven that 



the overlap integral is zero for this particular wall profile, (the fact that the moment 
eigenvalues are different in the broken and unbroken phases introduce a factor — 
which is forgotten in ref. fl^). Now we will present the general argument. 



The ultimate reason for the orthogonality of the eigenstates considered follows 
from the different physical processes to which they are related. The real processes are 
described by wave packets coming from the unbroken (broken) phase, bouncing on 
the wall, and generating reflected and transmitted packets. The solutions which are 
asymptotically plane waves have no real physical meaning, but they allow to build 
up localized wave packets. As we will see, the solutions named ■?/'^"'^ are associated 
to wave packets evolving from the unbroken phase to the broken phase and those 
named ip'^ describe the same but in the opposite direction, from broken to unbroken 
phase. Since these wave packets do not overlap in the past, by unitarity they will not 
overlap in all their evolution. Let us give a few more details. 

The X and y components of the wave packet are not modified by the wall, and 
will be ignored in what follows. We consider the following incoming wave packet, 
approaching the wall from the unbroken phase, 

which has the helicity, h, fixed. A^ is a normalization constant, the quantity r = 
z + Z — t — T has been introduced, where Z/d » 1 and T/d >> 1. d denotes the 
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spatial extension of the wave packet which is located at the time ~ —T around the 
position ~ —ZiZ > 0). d is introduced in such a way that p'^^^d » 1. It can be 
easily proven that at t ~ — T, 

P(r^;^,t) = N' f dhA{hJrWn%^), (81) 



where N' is a new normalization constant and n = kx, ky,jz, Ek, r. When the terms 
exponentially suppressed by e'^^'""'^^ are neglected as in the second line of (^OD, we 
obtain 



A(A;„pr) = e-(P^"^-^^)'^e^(^-*-^)'=^ (82) 

In other words, an incoming wave packet located far away from the wall in the 
unbroken phase, with group velocity pointing towards the wall, totally expands on 
the eigenstates called incoming. In |^ this general conclusion for any wall profile 
is stated for the particular case of the thin wall (step profile function). Analogous 
conclusions hold for the wave packets coming from the broken phase and for the 
eigenstates we called broken incoming (ipn)- Obviously the superposition of these 
two incoming wave packets located far from the wall, to the left and to the right, 
is zero at initial time — T. As they evolve in the time, for t = 0, they coincide 
in the same spatial region but their overlap integral keeps on being zero, since the 
Hamiltonian being hermitian the time evolution operator is unitary. It follows that 
the superposition of the two wave packets located in a certain region centered at 
z = 0, which extends a distance of the order of d, is zero. If the limit li — >■ cx) is 
taken now, the gaussian with the appropriate normalization constant gives a delta 
function in (|82D , and the overlap of the two wave packets becomes the overlap of ip^^'^ 
and V'^^. Obviously, d and Z go to oo so that Z/d ^ 1, our argument being valid 
in this limit. From the vanishing of that overlap we conclude to the orthogonality of 
these eigenf unctions. 
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1 



FOOTNOTES 

In these references the final expressions for R{E) and R{E) are equivalent to 
(^, although the notation is different. Moreover, in these references zq is taken 
to be 0. 

This is a reasonable assumption. By considering the appropriate dimensions we 



have that ^ = O[10GeV] (see for example [|T8|,0) and therefore a 



he 



O[10] in the electroweak phase transition energy scale. For this range of energy 
/' £ ^2 ^ ^ except for the top mass. 



^mo' 



3 



As known the function arctan(x) defined from R to (0, 27t) is bivalued. The 
right-hand of the second equation in (^) is well-defined by considering the sign 
of (F'^'^)(l) -|- F^''\—l)) sin a to determinate if the angle belongs to (0, tt) or to 
(tt, 27r). If it is zero, {F^^\l) — F^^\—l)) cosa must be considered in order to 
determine if the angle is or tt. 



5 



6 



In ref. [Q we use analytic profiles which behave as e'"~^ " , 6"*"+^ " , for 
X << —1, X >> 1, respectively. The results which will be derived in what 
follows can be, in the same way, obtained for these profiles. 

The position of the closest pole to the real axis, i/j, gives us the smallest value 
of A we can consider in order to be able to neglect the contribution of Rn-c- 
We have A ^ ^^Vj^ ^§^yj- ^^ ^^ ^^^^ ^^ ^^^^ ^^^ term Rn-c is associated 
to the fact of considering the evolution in the region (— 5^„, 5yj) and introducing 
the technical assumption mentioned above. 

In order to compare results the relationship a = 1/a, where a is the energy 



parameter used in ref. [IS], must be taken into account. 
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